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Abstract 

In this work, we study the evolution of the susceptible individuals during the spread of an 
epidemic modeled by the susceptible-infected-recovered (SIR) process spreading on the top of 
complex networks. Using an edge-based compartmental approach and percolation tools, we find 
that a time-dependent quantity 3>s(t), namely, the probability that a given neighbor of a node 
is susceptible at time t, is the control parameter of a node void percolation process involving 
those nodes on the network not-reached by the disease. We show that there exists a critical 
time t c above which the giant susceptible component is destroyed. As a consequence, in order to 
preserve a macroscopic connected fraction of the network composed by healthy individuals which 
guarantee its functionality, any mitigation strategy should be implemented before this critical 
time t c . Our theoretical results are confirmed by extensive simulations of the SIR process. 



Electronic address: ldvaldes[at]mdp.edu.ar 



I. INTRODUCTION 



The study of epidemic spreading has been one of the most successful applications on 
networks science. Recent outbreaks of new influenza strains like the H1N1 |l| and the 



H5N5 flu or the Severe Acute Respiratory Syndrome (SARS) |2j, which are characterized 
by a high rate of mortality and/or fast propagation velocity, motivate the development 
of epidemic models that capture the main features of the spread of those diseases. In 
particular, mathematical tools applied to model epidemics are very important since they 
allow to understand how a disease impact on the society, helping to develop new policies 
to slow down its spreading. 

One of the simplest models that reproduce seasonal diseases, such as influenza, is the 
susceptible-infected- recovered (SIR) model 3|, |4[, which has been the subject of extensive 
theoretical and numerical research on complex networks In the SIR model the indi- 
viduals can be in one of three states, susceptible, infected or recovered. In its discrete 
formulation p-L7j, at each time step, infected individuals infect their susceptible neighbors 
with probability (3 and recover at a fixed time t r since they were infected, called recovery 
time. According to these rules, the disease spreads on the contact network until it reaches 
the steady state where there are only susceptible and recovered individuals. It was found 
that the steady state of the SIR model can be mapped into a link percolation problem 
which provides a theoretical framework to study this process JgI-IsUiOI. It is known that the 
size of the infection, defined as the fraction of recovered individuals at the steady state, is 
governed by the effective probability of infection or transmissibility T of the disease which 
depends on f3 and t r . In the SIR model, the size of the infection is the order parameter 
of a second order phase transition with a critical threshold transmissibility T c . Below T c 
the disease is an outbreak, where the infection reaches a small fraction of the population 
while above T c an epidemic develops exactly as in a link percolation process J BQ. 
In uncorrelated infinite networks this threshold is given by T c = 1/0 - 1) |6|, Hill, where 
k = (k 2 ) I (k) is the branching factor of the network, and (k) and (k 2 ) are the first and 
the second moment, respectively, of the degree distribution P(k). Here, k is the degree 
or number of links that a node can have with k m i n < k < k max . For Erdos-Renyi net- 
works (ER), the degree distribution is P(k) = (k) k /k\ and the threshold is found at 
T c = 1/ (k) . However, most of the real networks have a heterogeneous degree distribution 
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that is better represented by a pure Scale-Free network (SF) with P(k) ~ k~ x , where A 
measures the broadness of the distribution. In the thermodynamic limit, for SF networks 
with 2 < A < 3, (A; 2 ) — > oo and as a consequence, the critical transmissibility T c — > 
which means that the epidemic spreads for any value of T . However, due to finite 

size effects, real networks have finite critical transmisibilities. 

In a recent paper, using a generating function formalism, Newman 12[ showed that 



at the steady state of the SIR model there exists a second threshold T* above which 
the residual network composed by the biggest giant susceptible cluster that remains after 
a first propagation, is destroyed. From an epidemiological point of view, this implies 
that if a disease spreads for a second time on the residual network, it cannot become 



13] showed that T* is an important 



an epidemic. On the other hand, Valdez et al. 
parameter to determine the efficiency of a mitigation or control strategy, because any 
strategy that decrease the transmissibility below T*, can protect a large and connected 
cluster of susceptible individuals. Using a percolation framework, they explained the lost 
of the susceptible giant cluster as a not-random node percolation process, that they called 
node void percolation, in which a susceptible individual corresponds to a void node in link 
percolation. 

Even though percolation theory was very useful to describe the steady state of the 
SIR model on complex networks, it is still very challenging to explain the dynamics 
of the model to develop intervention strategies before the epidemic spreads to a large 
fraction of the population. To describe the dynamics of epidemic spreading on networks, 
recently some researchers developed differential rate equations for the SIR model that 
take into account the network topology. Lindquist et al. [3| introduced an "effective 
degree" approach through a large system of ordinary differential equations. Under this 
approach, the nodes and their neighbors are categorized by their disease state (susceptible, 
infected, recovered) and each differential equation compute the evolution of the fraction of 
susceptible or infected nodes with a number i and s of infected and susceptible neighbors, 
respectively, with < % < k max and < s < k max . As a result, a system with O (A^ aa .) 
equations needs to be solved. This approach represents accurately the evolution of the 
number of infected individuals, but at a high computational cost. On the other hand, 
Miller 15 1 and Miller et al. flf] . \v\ proposed an ingenious approach to describe the 



evolution of a SIR process with rates by means of an edge-based compartmental model 
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(EBCM) |l5l . Il6] which has the advantage to describe the dynamical spreading of an 
epidemic with only a few equations. With these equations, the authors found accurate 
results for the evolution of the number of infected individuals for static and dynamic 
evolutive topologies like "edge swapping" and "dormant contacts" for transmissibilities 
above the critical threshold |16| . 

While most of the literature is focused on studying the evolution of the fraction of in- 
fected or susceptible individuals, it has not yet been investigated how the epidemic spread 
affects the evolution of the network composed by the susceptible individuals. Understand- 
ing this problem is important because the network composed by the healthy individuals 
is the network that sustains the functionality of a society, e.g. the economy of a region. 
In this paper we present a novel idea for the SIR model, based on a dynamical study of 
the network composed by susceptible individuals. We show that the temporal decreas- 
ing of the size of the giant susceptible cluster can be described as a dynamic void node 
percolation process with an instantaneous void control parameter. We find that there 
exists a critical time t c above which the giant susceptible component overcomes a tem- 
poral second order phase transition with mean field exponents. The paper is organized 
as following: in Methods and Results we present the theoretical framework to derive the 
evolution equations. Then we study the evolution of the giant susceptible cluster and its 
temporal critical behavior. Finally we present our conclusions. 

II. METHODS AND RESULTS 
Theoretical framework 

The evolution equations of the dynamic SIR model provide the basis for analyzing 
theoretically novel magnitudes that could be useful for epidemiologists and authorities to 
plan policies to stop a disease before an epidemic develops. In the SIR model, initially, 
all the nodes are susceptible except for one node randomly infected, that represents the 
index case from which the disease spreads. The infected individual transmits the disease to 
susceptible neighbors with probability /3 and recovers t r time units since he was infected. 
For the SIR with fixed recovery time, the transmissibility is given by T(f3,t r ) = T = 



1 - (1 (l3|. 
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In order to study the evolution of the states of the individuals in the SIR with fixed 



recovery time, we use the edge-based compartmental model (EBCM) 15Nl7|. The EBCM 



is based on a generating function formalism widely implemented in branching and per- 
colation process on complex networks j^l, [isL 20 1 . For a branching process that spreads 
on uncorrelated networks, such as the tree of infected individuals, two generating func- 
tions that contain the information of the topology of these networks are defined. The 
first one is the generating function of the node degree distribution P(k) which is given by 
Gq(x) = P(k)x h . The second one is the generating function of the degree distribution 
of the first neighbors of a node, also called excess degree distribution Pi(k) = kP(k)/(k), 
given by G\(x) = kP(k)/ (k)x k ~ l . Here, P\(k) is the probability to reach a neighbor 
of a node, following a link. It is straightforward that the mean connectivity of the nodes 
is (k) = G' (l). 

Denoting the fraction of susceptible, infected and recovered individuals at time t by 
S(t), I(t) and R(t), respectively, the EBCM approach describes the evolution of the 
probability that a node (which we call root node) is susceptible. In order to compute this 
probability, an edge is randomly chosen and a direction is given, in which the node in the 
target of the arrow is the root, and the base is its neighbor. Disallowing that the root 
infects the neighbor, 8(t) = 9 t is the probability that the neighbor does not transmit the 
disease to the root, with 9 t given by 

6 t = $ s (t) + $ I (t) + $ R (t), (1) 

where $s(t), <&#(£) and are the probabilities that the neighbor is susceptible, recov- 

ered, or infected but has not transmitted yet the disease to the root. The probability that 
a root node with connectivity k is susceptible is therefore Q\ and the fraction of susceptible 
nodes is S(t) = Ylk^i^^t = Go{&t)- This approach simplifies the calculations, reducing 
the problem to finding an evolution equation for 9 t , from where the evolution of S(t), 
R(t) and I(t) is derived. Thus, using the EBCM approach adapted to SIR with fixed t r 
(see Appendix Sec. I), the evolutions of 6 t , &s{t) an d &i{t) are given by the deterministic 
equations 
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A6 t = (2) 
A$ s (t) = G 1 (9 t+1 )-G 1 (9 t ), (3) 
A$j(t) = - A$ s (t) + (1 - T)A$ s (t - t r ), (4) 

where A is the discrete change of the variables between times t and t+1. Eq. ([2]) represents 
the decrease of 9 t when a infected neighbor transmits the disease. Eq. represents the 
decrease of $s(t) when a susceptible neighbor is infected (notice that A$s(t) < 0). This 
term contributes to an increase of in Eq. (j3J) where the first term represents the 

decrease of $/(£) when the links transmit the disease, the second term corresponds to the 
term of Eq. ([2D mentioned above and the third term represents the decrease of due 
to the recovery of infected individuals. 

From the above equations, the evolution of the fraction of infected individuals can be 
computed as 

AI(t) = -AS(t) + AS(t-t r ), (5) 

where the first term represents the fraction of new infected individuals (see Appendix 
Sec. I). The second term represents the recovery of infected individuals that have been 
infected t r time units ago. 

These difference equations correctly describe the evolution of S(t), I(t) and R(t) above 
the criticallity for all values of t r and /3 (see Appendix Sec. I). In the next section, we will 
show that combining this approach and dynamic percolation, we can describe the time- 
dependent evolution of the susceptible individuals in the SIR model as a dynamic void 
node percolation process for any value of t r . 



A. Temporal percolation of susceptible individuals 



In Ref. 13| it was found that the process under which the susceptible clusters size 
decrease can be explained with node void percolation defined below that as we will show 
can be related with the dynamic SIR process. 

In the steady state of the SIR model an epidemic cluster is equivalent to a Leath 
growth process J21I Q| with a link occupancy probability T. The Leath process on 
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complex networks generates a single cluster that represents the infection tree for a given 
value of the transmission probability T. Denoting by f n (T) the probability that a cluster 
reaches the nth generation following a link, the probability foo{T) that a link leads to a 
giant component {n — > oo) is given by 
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22] 

/oo(T) = l-^^[l-T /oo (T)f- 1 , (6) 
k=i W 

where /oo(T) is the solution of 

f OD (T) = l-G 1 [l-T U(T)}. (7) 

As the "infectious" cluster grows from a root, generation by generation, the sizes of the 
void clusters, i.e. the nodes not reached by the disease, are reduced as in a node dilution 
process, since when a link is traversed a void cluster loses a node and all its edges. As a 
consequence, for large generations foo{T) can also be interpreted as the probability that 
a void cluster loses a node. However, in this kind of percolation process the void nodes 
are not killed at random, instead they are removed following a link. We call this type of 
percolation "node void percolation". If we denote by 1 — V s the probability that a void 
node is removed due to the occupancy of a link, at the steady state the following relation 
holds 

l-^ = /co(T). (8) 

Then V s is the probability that a void node is not removed due to the fact that the 
link has not been traversed. Thus, V s is equivalent to $s(£ — > oo) because the void 
nodes correspond to the susceptible individuals in the steady state. As in any percolation 
process, there is a critical probability V£ at which the void network undergoes a second 
order phase transition. Above V£ a giant void component exist while at and below V£ 
void nodes belong only to finite components. In epidemic terms, this means that at V* 
only finite susceptible clusters can be reached. As a consequence, the fraction of links T* 
needed to reach this point fulfills {3| 

K s = Woo(T*). (9) 

Therefore, from Eqs. (J2j) and we obtain 

V^G^l-Til-Vi)], (10) 



where T* is the solution of Eq. (TTU]) . This result shows that at the steady state, for T >T 



we have V s < V c s and therefore the size of the giant susceptible cluster Si — > 



13|. Even 



though static percolation is a useful tool to analyze the final size of the giant component 



of susceptible individuals 



12(| , it is very important to know the evolution of S\(t), since it 



can be used as a criteria to begin or to increase an intervention to protect a large fraction 



of the susceptible population jl3|. As we will show below, Si(t) can be fully related with 
a node void percolation process at every instant t. 

In order to describe the evolution of the size of the giant susceptible cluster, we define 
bj t as the probability that a neighbor of a root not connected to the giant susceptible 
cluster has not yet transmitted the disease to the root at time t. This is possible if the 
neighbor of the root node is infected but has not yet transmitted the disease, recovered 
or susceptible but not connected to the giant susceptible cluster, with probabilities $/(£), 
$_r(£) and Gi(u t ) respectively. Similarly to 9 t (see Eq. ([I])), these probabilities satisfy the 
relation 

* R {t)+$ I (t)+G 1 (L> t )=U t , (11) 

where Gi(u t ) is the generating function of the neighbor of a root not connected to the 
giant susceptible cluster. From Eq. ((U), $ 7 (t) + § R (t) = 9 t - $ 5 (t) = 9 t - Gi(9 t ). Then 
Eq. (ITT]) can be rewritten as, 

ut-Gxiut) = 9 t -G 1 (9 t ), (12) 

and the evolution of Si(t) is given by 

Si(t) = G (9 t )-G (cu t ), (13) 

where Go(9t) is the total fraction of susceptible individuals and Go(ut) is the fraction of 
individuals belonging to finite susceptible clusters at time t. Notice that the dynamical 



Eqs. (fl2|) and ( JT3|) are a time-dependent versions of the ones derived in Ref. [12j for the 
steady state (t — > oo) of the SIR model. This suggests that the evolution of the giant 
susceptible or percolating void cluster can be thought as a temporal percolation process. 
Thus, the magnitudes derived for the static percolation of the susceptible individuals have 
a dynamical counterpart. As a result, V s and &s{t), are equivalent not only at the steady 
state, but also at every instant of time. In order to show the equivalence, in Fig.[T]we show 



in the same plot Si(t) as a function of $s(t), obtained from Eqs. ©-([2]) and (fl2"|) -( 1T3"|) . 
and the steady state S\(t — > oo) as a function of V s 12| for ER and SF networks with 
the same (k) and iV for T = 0.76 > T*. 




FIG. 1: Equivalence between $s(t) and V s . S\(t — > oo) as a function of V s (x) obtained 



in Refs. 
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and Si(t) as a function of (solid line) obtained from Eqs. ©-(ED 

and (|T2|) -( JT3|) with N = 10 5 and mean connectivity 4.07 in the giant component for (A) a 
ER network with (k) = 4 and (B) SF network with A = 2.63, k m i n = 2 and (k) = 4.07. In 
the insets we show Si(t) as a function of $s(t) from the simulations (symbols) and from 
Eqs. ©-([2D and (fl2>(jI5]) (solid line) for t r = 1 (□) and t r = 20 (Q). (Color online). 



As we can see, the static curve Si(t — > oo) as a function of V s is the same as Si(t) 
as a function of $s(i) and they coincide with the simulations for different values of t r 
which shows the equivalence between V s and &s{t) & t every instant of time and not only 
at the steady state (for details of the simulations see Appendix Sec. I). Thus our process 
can be explained by a dynamic percolation with an instantaneous void transmissibility 
V s = ® s (t). 

With our theoretical formulation, we will show that there is a critical time t c at which 
the giant susceptible cluster disappears that correspond to the time at which $,s(t c ) = V/. 
In order to prove this, notice that according to Eq. (1T21) . 9 t and oj t can be thought as two 
points with the same image of the function x—G\ (x) . Solving this equation for the variable 
uj t above T*, two solutions are found since the curve x — Gi(x) is a concave function for 
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x > as can be seen in Fig. [2j One of the solutions is the trivial one, for which S\(t) = 0, 
that corresponds to the maximum of the function x — G(x) at 6 tc = u tc = u c . Then the 
giant susceptible cluster is destroyed at the point uj c which fulfills 



[x-Gx(x)] 



0. 



(14) 



then, 



G[) (1). 



(15) 



Thus when Eq. (TT4"j) is satisfied, the giant susceptible cluster disappears and $s(i c ) = 
V' = GAu c = 9 C ) = G x \(G\Y X (i) , i.e. 



(G , 1 r(i) j 

®s(tc) = G 1 



For ER networks it is straightforward to show that §s{t c ) = 



(16) 



I 

X 




X 



FIG. 2: Schematic of the behavior of Eq. fTT2|) for T > T*. From the initial condition 
#o = 6(t = 0) = 1, 9 t and uj t , satisfies Eq. (|T2|) . For 9 t ^ oJt we have two solutions 
that correspond to S\(t) > 0. When 9 t reaches the maximum of the function x — Gi(x), 
6 C = u c , the giant susceptible component is destroyed. The dashed lines are used as a 
guide to show the possible solutions of Eq. (I12p . 
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In Fig. [3] we plot the time evolution of the fraction of susceptible individuals Si(t) in 
the susceptible giant component as a function of t for ER and SF networks obtained from 
the theory and the simulations, for a transmissibility T above T*. 




FIG. 3: Time evolution of Si(t) for t r = 20 and /3 = 0.07 (T = 0.76) and mean connectivity 
4.07 in the giant component for (A) a ER network with (k) = 4 (T* = 0.46) and (B) a 
SF networks with A = 2.63, minimal connectivity k m i n = 2 and (k) = 4.07 (T* = 0.38). 
The symbols correspond to the simulations with the time shifted to t — when 1% of 
the individuals are infected, and the solid lines correspond to the theoretical solutions 
Si(t) (blue solid line) of Eqs. (|T2|) - (fT3|) . In the insets we show the size of the second 
biggest susceptible cluster S^t) ( re d solid line) and the evolution of $s(t) (black solid 
line) obtained from simulations. The value of $s(£ c ) = V* (dashed line) was obtained 
from Eq. (fT6|) . ^(t) has been amplified by a factor of 50 in order to show it on the 
same scale as the rest of the curves. The simulations are averaged over 1000 network 
realizations with N = 10 5 . (Color online). 

As shown in Fig. El there is an excellent agreement between the theoretical curve Si (t) , 
obtained from Eqs. (IT2"|) and (I13p . and the simulations which validate that percolation 
tools can be used to describe the time dependence of the susceptible individuals in the 
SIR process for T > T*. On the other hand, in the figure we can see that for T > T*, the 
giant susceptible cluster Si(t) is destroyed at t = t c which occurs exactly at $s(£ c ) = V* 
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(see the insets of Fig. [3]). Our results show that can be used to determine whether a 

giant susceptible cluster exists at a given time. In turn, in the insets of Fig. Owe can see 
that the size of the second susceptible cluster S%(t) has a sharp peak around t c , indicating 
that, as in static percolation, the susceptible individuals overcome a second order phase 
transition. However, this transition is not given by a random node percolation process. 
As the disease spreads through the links, the susceptible individuals are removed with 
probability proportional to kP(k), i.e., the susceptible network loses the higher degree 
nodes first. For this reason, the disease spreading induces a second order phase transition 
in the susceptible network with mean field exponents at t c (see discussion in the Appendix 
Sec.II). 

An important implication of our results is that, it can be used by the health authorities 
to implement intervention strategies before the critical time t c is reached. This will allow 
to protect a macroscopic fraction of the network composed by healthy interconnected 
individuals which preserve all the topological properties characteristic of social contact 
networks and their functionality. 

Conclusions 

In this paper we introduce a temporal dynamic percolation to characterize the evolu- 
tion of the susceptible individuals in a SIR model. We show using an edge-based com- 
partmental approach and percolation tools that as the disease spreads the evolution of 
the susceptible network can be explained as a temporal node void percolation that can be 
mapped instantaneously into static percolation. We show that for transmissibilities above 
T*, there exist a critical time above which the giant susceptible cluster is destroyed and 
the susceptible network overcomes a second order transition with mean field exponents. 
All our theoretical results are in excellent agreement with the simulations. Our findings 
are very interesting from an epidemiological point of view since the existence of a thresh- 
old time implies that when a very virulent disease reaches a small number of susceptible 
individuals, the authorities have only a limited time to intervene, in order to protect a 
big community (susceptible giant component) that has not been already reached by the 
epidemic, and to preserve the topological features of SF networks. Our finding on the 
susceptible network could be extended to other epidemics dynamics allowing to obtain a 
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better description of the effect of diseases spreading on social and technological networks. 
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Appendix A: 


The EBCM approach 




TABLE SI: Definitions 


Variable/Parameter Definition 


P 


Infection rate or infection probability. 


7 


Recovery rate. 


U 


Recovery time. 


T 


Transmissibility. 


o, 


Probability that a neighbor of a root node has not transmitted yet the 




disease to the root node at time t. 




Probability that a neighbor of a root node is susceptible at time t. 




Probability that an infected neighbor of a root node has not transmitted 




the disease to the root node at time t. 




Probability that a neighbor is recovered at time t without having trans- 




mitted the disease to the root node. 



In the EBCM approach, Q t is the probability that a root node has not being infected by 
a neighbor at time t. This is possible if the neighbor is susceptible, recovered, or infected 
but has not transmitted the disease yet to the root, which happens with probabilities 
$s(t), and $/(*), respectively. Then, 9 t = $ s (£) + $/(t) + $ R (t). The probability 

that a root node of connectivity k is susceptible is 9 k and the fraction of susceptible 
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nodes is S(t) = P(k)6^ = Go(6 t ). On the other hand, a neighbor is susceptible 
with pro bability $s(t) = Gi(8 t ). Then in the SIR model with infection and recovery 
rates 15l Il6|. the probabilities $/(t), $s(£) an d 6t evolve as, 



6 = (SI) 
$ 5 = -PG'^O)®!, (S2) 
$/ = -/?$/ + - 7 $ 7 , (S3) 

where /3 and 7 are the infection and recovered rates. Eq. (1S1I) represents the decrease of 
# when an infected neighbor transmits the disease. The Eq. (1S2|) represents the decrease 
of $5 when a susceptible neighbor is infected, which is proportional to G'^O), i.e., the 
mean connectivity of the susceptible first neighbors or the excess degree of the susceptible 
individuals, because when a susceptible individual is infected, all its links except the one 
used to infected it, can transmit the disease. This term contributes to an increase of 
in Eq. (1S3I) . In Eq. (jS3l) on the r.h.s, the first term represents the decrease of when the 
links transmit the disease, the second term corresponds to the term of Eq. (1S2I) mentioned 
above and the third term represents the decrease of due to the recovery of infected 
individuals. 

To obtain the evolution of I(t), we use the fact that, 

I + S + R = 0. (S4) 

As R = 7 / and S = d (G (6)) /dt = G' (6)9 = -/3$/G' o (0), the evolution of the fraction 
of infected individuals is given by 

i = pCf^j—yl, (S5) 

where the first term represents the decrease of S which is proportional to 0, the mean 
connectivity of susceptible individuals G' (9) and the probability that an outgoing edge 
from a root is connected with an infected node that has not transmitted the disease to 
the root at time t. The second term corresponds to the recovery of infected individuals 
at a rate 7. 

We reformulate the EBCM approach process with discrete time steps, for a fixed re- 
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covery time t r . It is straightforward that Eq. (jS 1HS5[) can be written as, 

A9 t = -/3^j(t), (S6) 
A$ s (t) = G l (e t+1 )-G 1 (e t ), (S7) 
A$j(t) = r(t) - A$ s (t) + (1 - T)A$ s (t - t r ), (S8) 

where 1 — T = (1 — /3)* r denotes the probability that an infected individual has not 
transmitted the disease to a susceptible individual during t r time units since he was 
infected. Finally the evolution of the fraction of infected individuals is given by 

AI{t) = -AS(t) + AS(t-t r ), (S9) 

where —AS(t) = — (Go(#t+i) — Go(0t)) represents the fraction of new infected individuals 
and the second term represents the recovery of infected individuals that have been infected 
t r time units ago. 

For the simulations we infect only one individual in the giant component of the network 
and at each time step all the infected individuals infect their susceptible network with 
probability /3 and recover at a fixed time t r since they were infected. We select only the 
runs in which the size of the epidemic has reached a macroscopic fraction of individuals in 
the steady state 5|, |6|, |9| because the deterministic equations are only valid for epidemics 
above the critical threshold T c . We performed all the simulations using synchronized or 
simultaneous updates at each time step. 

In Fig. li?Tj we plot the time evolution of the fraction of infected nodes I(t) for ER 
and SF networks obtaining by the EBCM approach Eqs. (1S6HS9|) and the simulation. For 
the simulations we shifted t = to the instant when the disease has reached 1% of the 
individuals. We choose this reference time, as the time when the disease has reached a 
size enough to growth deter ministically. This choice compensates the time dispersion of 
each trial around the theoretical solution due to stochastic effects at the early stages of 
the process when the number of infected nodes is small 16| (see the insets of Fig. IS 1 [) . As 



shown in Figs. ISTR-B. each trial simulation has the same shape as the theoretical solution 
which shows that the EBCM approach and the simulations are in excellent agreement. 
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FIG. SI: I(t) for epidemics with t r = 20 and f3 = 0.04 (T = 0.55) on networks with 
mean connectivity 4.07 in the giant component, for a ER network with (k) = 4 (A) and 
a SF with A = 2.63, minimal connectivity fc min = 2 and (k) = 4.07 (B). The symbols 
correspond to an average of one hundred different network realizations with N = 10 5 
nodes and the solid black curve is the numerical solution of Eq. (1S9I) shifting the curves 
to t = when 1% of the individuals are infected. The insets show the individual 100 
network realizations (solid gray lines) and the numerical solution of Eq. flS9j) (solid black 
line) without the temporal shift transformation. 



Appendix B: Node void percolation in the time domain 

In node void percolation, as a link is traversed, void node is removed. The void 
nodes are removed with probability proportional to kP(k). As the susceptible nodes can 
be mapped into void node percolation, the susceptible network loses their higher degree 
nodes first as in an intentional attack. As a consequence, the resulting susceptible network 



is more homogeneous than t 



re original. Thus, mean field exponents of a second order 



percolating phase transition 23| are expected. In order to show the effect of the disease 
spreading on the highest degree nodes, in Fig. [S2] we plot for a SF network the effective 
degree distribution of the susceptible nodes obtained from the simulations, in which a 
susceptible node has degree k when it has k susceptible neighbors. 
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FIG. S2: Simulation results of the degree distribution of susceptible nodes for a SF network 
with A = 2.63, k m i n = 2 and (k) = 4.07 at different times: at the beginning of the 
spreading (black solid line), when the disease has reached 10% of individuals (blue solid 
line), 25% of individuals (violet solid line) and 50% of individuals corresponding to t c (red 
solid line). (Color online). 

As shown in Fig. IS2| as the disease spreads, the effective degree distribution loses the 
heavy tail. As a result of this process, the susceptible clusters becomes more sparse and 
at the critical time t c the topology of the susceptible clusters change drastically since the 
susceptible individuals lose all the hubs and P(k) has an exponential tail. For percolation 
in mean field it is known that at the criticality the finite cluster size distribution n$ ~ s~ T 
with t = 2.5 and Si [$ s (^)] ~ &s{t) — & s {t c )- in Fig. [S3]we plot the simulations results of 
the finite size distribution of the susceptible nodes n s at t = t c . 
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FIG. S3: Log-log of the cluster size distribution n s of finite susceptible clusters (O) & t t c 
for t r = 20 and (3 = 0.07 (T=0.76) in a ER network with (k) = 4 for t c = 30(A) and a SF 
with A = 2.63, minimal connectivity k m i n = 2 and (k) = 4.07 for t c = 11. (B). The dashed 
line corresponds to a power law fitting, from where we obtain an exponent r w 2.5. Our 
simulations were averaged over 10000 network realizations with N = 10 5 . 
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FIG. S4: S\ as a function of 3> s (t) obtained from simulations (O) an d from the analytical 
approach (solid line) in a ER network with (k) = 4 (A) and a SF with A = 2.63, minimal 
connectivity k m i n = 2 and (k) = 4.07 (B). In the inset we plot Si as a function of 
the distance of $ s (£) to the criticality $ s (£ c ) = V'c, in log-log scale. The dashed line 
corresponds to a power law fitting from where we obtain slope ~ 1. Our simulations were 
averaged over 1000 network realizations with N = 10 5 . 



We can see that at t c , n s (t c ) behaves as a power law with exponent r ^ 2.5 which 
corresponds to the mean field value, independently of the initial degree distribution of 



the network 



24j . Similarly, in Fig. ED we plot Si(t) as a function of $s(t) obtained from 



the simulations and the theoretical approach. We compute $s(t) from the simulations as 
the square root of the fraction of edges connecting two susceptible nodes 25]. We can see 
that Si(t) behaves as a power law with exponent one with the distance to the criticality 
*&s{tc), which also corresponds to the mean field value (see Insets of Fig. 154]) . Since two 
critical exponents are sufficient to determine the universality class, the results showed 
above indicate that in a node void percolation process the susceptible network belongs 
to the same universality class of mean field percolation and confirms quantitatively the 
homogenization of the susceptible network during a SIR epidemic spreading. 

Finally, in Fig. [55l we plot the critical time t c , computed from the simulations at 
*&s{t) = Vc-i as a function of T for different values of t r . We can see that for the same 
transmissibility T, when t r increase, the time to intervene grows since /3 decrease and 
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thus the disease spreading is retarded. In turn, when the transmissibility T reaches 
T* from above, the critical time t c grows very fast. This phenomenon is analogous to 
other second order phase transitions in physics like the relaxation time near the Curie 
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27j, and indicates that that 



temperature, which are called "critical slowing down'' 
once the transmissibility increases slightly above T*, the time needed to destroy the giant 
susceptible cluster decreases very fast. 
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FIG. S5: t c as a function of T for /3 = 0.07 and t r = 20 (O), U = 10 (□), t r = 5 
(o) and mean connectivity 4.07 in the giant component in a ER network with (k) = 4 
(T* = 0.46) (A) and in a SF with A = 2.63, minimal connectivity k min = 2 and (k) = 4.07 
(T* = 0.38)(B). The dashed line represents the value of T*. The critical time t c is 
measured using t = when 1% of individuals are infected. The dotted lines are used as 
a guide to the eyes. 
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